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CN . Abstract 

Relativistic quantum transport theory has begun to play an important role in the space-time 
description of matter under extreme conditions of high energy density in out-of-equilibrium situ- 
ations. The following introductory lectures on some of its basic concepts and methods comprise 
the sections: 1. Introduction; 2. Aims of transport theory (classical); 3. Quantum mechanical 
distribution functions - the density matrix and the Wigner function; 4. Transport theory for quan- 
tum fields; 5. Particle production by classical fields: 6. Fluid dynamics of relativistic quantum dust. 

(N 
O 

P-i. 1 Introduction 

The study of the behavior of matter under more and more extreme conditions has a long tradition 
motivated by the quest for understanding the forces among its constituents on smaller and smaller 
scales. Not only the attempts to understand such spectacular phenomena as the stellar supernova 
^ explosions or theories of even the primordial stages of cosmological evolution, but also the ever in- 

creasing collision energy of high-energy particle accelerators and the heavy-ion programs at CERN 
and RHIC, in particular, witness the most recent stages of this scientific development. 

In the latter context, often the complicated space-time dependence of what is really a quantum 
many-body system or what are highly dynamical interacting quantum fields is described in terms of 
a perfect fluid model. Since the seminal work by Fermi and Landau this approach has been applied 
successfully, in order to study global features, such as multiplicity distributions and apparently thermal 
transverse momentum spectra of produced particles, in high-energy collisions of strongly interacting 
matter ||, ||, || Similarly, the hydrodynamic approximation is often invoked in astrophysical 
applications and cosmological studies of the early universe || 

The limitations of and likely necessary corrections to the fluid picture, however, have rarely been 
explored in the microscopic or high energy density domain. Difficulties reside in the derivation of 
consistent transport equations and in the amount of computational work required to find realistic 
solutions; see Refs. §| 0, ||, for example, for a review and recent progress concerning selfinteracting 
scalar particles and the quark-gluon plasma, respectively. More understanding of related hydrody- 
namic behavior, if any, seems highly desirable. 

For example, it has recently been shown that a free scalar field indeed behaves like a perfect 
fluid in the semiclassical (WKB) regime |]]. More generally, the mechanisms of quantum decoherence 
and thermalization in such systems which can be described hydrodynamically, i.e. the emergence of 
classical deterministic evolution from an underlying quantum field theory, are of fundamental interest 
[1, |1 0, 11. 
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Having said this, it becomes obvious that not only particular applications of relativistic quantum 
transport theory motivated by experiments or observations - such as a suitable quark-gluon plasma 
transport theory to be applied to the phenomenology of high-energy heavy-ion collisions, or a transport 
theory for the electro-weak interactions of the intense neutrino flux from a supernova core with its 
electron-positron plasma sphere - are of interest, but that many interesting conceptual problems can 
be found in this field. It is the aim of the present introduction to describe some of its basic concepts 
and methods. 

It seems worth while to emphasize here that transport theory by its very nature aims to describe 
highly dynamical systems where the time dependence of the phenomena to be studied cannot be ne- 
glected. Therefore, one necessarily has to go beyond (thermal) equilibrium field theory, for example. A 
partial exception consists in linear response theory, where standard field theory methods are employed 
to calculate the short-time response of the system to necessarily small perturbations. 

The plan of these lectures is as follows. In Section 2 we present the motivation for transport theory 
by taking a cursory look at classical relativistic transport theory and its relation to relativistic hydro- 
dynamics. In Section 3 we introduce quantum mechanical distribution functions. Especially, the need 
for the density matrix formalism is reviewed in basic terms, and the Wigner function is introduced. 
In Section 4, as an example, we develop the transport theory for the particular model of interacting 
quantum fields with a global 0(4) symmetry, i.e. the linear sigma model, in the Hartree approxima- 
tion. In Section 5 particle production by classical fields is described, presumably an important effect 
during early stages of heavy-ion collisions, solving fermion quantum transport equations perturba- 
tively. Finally, in Section 6, we investigate the fluid dynamical behavior of relativistic quantum dust, 
solving the free quantum transport equations for arbitrary initial conditions exactly. 

2 Aims of transport theory 

Concerning the historical development as well as the systems under study, the subject matter of 
transport theory is most frequently associated with nonequilibrium plasmas of all sorts: 

• The quark-gluon plasma with QCD interactions formed shortly after the Big Bang initiating 
the observable Universe and possibly recreated during high-energy collisions, in particular with 
heavy nuclear projectiles/targets (R+Q). 

• The vve + e~ plasma with electro-weak interactions created during supernova explosions between 
the proto-neutron star core and the leftover outer layers of the collapsing star (R) . 

• H, He, . . . fusion plasmas in burning stars. 

• Electrodynamic plasma phenomena in the Earth's ionosphere leading to polar lights and thun- 
derstorms with lightnings. 

• Discharge plasmas used in neon lighting and plasma welding. 

• The e~ plasma in metals or semiconductors, the study of which has been advanced in solid state 
physics with particular attention to quantum effects (Q). 

Here we marked by R and/or Q the systems where relativistic and/or quantum effects are known 
or expected to play an essential role. - Abstracting from these examples, we notice two common 
qualitative features among these systems: 

• microscopic distance scales (average interparticle distance n -1 / 3 , mean free path X m f p , etc.) 
pa homogeneity scale L; 
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• microscopic time scales (average lifetime 77, relaxation time(s) T r , etc.) 
~ hydrodynamic time scale Th = L/c s , where c s denotes the sound velocity. 

Clearly, there will be exceptions to these qualitative statements and more precise characterizations 
of the plasma state can be given in the individual cases. However, this may suffice here, and we 
embark on the more formal description in the following. - We refer the reader to the monograph 
on relativistic kinetic theory of Ref. [14], which presents an excellent detailed exposition of the more 
traditional material in this context. 



2.1 Classical phase space description of many-body dynamics 

For a classical many-particle system, consider the probability to find a particle in the 8-dimensional 
phase space volume element at the (four-vector) position x with (four-vector) momentum p, 

dP(x,p) = f(x,p)d 4 xd 4: p , (1) 

where f(x,p) denotes the corresponding Lorentz scalar phase space density. Generally, in a classical 
system, the four-momentum is on-shell, such that not all components of p are independent. In par- 
ticular, we assume that the constraint expressing the energy in terms of the three-momentum and 
particle mass, p° = + y / p 2 + m 2 , is incorporated in /. - Except when explicitly stated, our units are 
such that h = c = ks = 1, and we use the Minkowski metric g^ v = diag(l, —1, —1, —1). 

A remark is in order here. Clearly, in order to learn not only about the behavior of a typical 
particle but also about its correlations with others, one principally should study the one-body density 
/ along with a two-body distribution f(xi, X2,pi,P2), three-body distribution . . . , etc. This is quite 
complicated in general and we restrict our attention to the one-body density here. Sometimes it is 
useful to visualize this function alternatively as describing a collection of (test) particles or a single 
particle with an ensemble of initial conditions. 

Now, how does the one-body density / evolve, e.g. from one time-like hypersurface to another? We 
recall two ingredients of Liouville's theorem, which will provide the basic tool to answer this question 
1 15]: i) Consider a phase-space volume element which is defined by 'tracer' particles forming its surface; 



then, due the uniqueness of the Newtonian motion, or its relativistic generalization, the number of 
particles inside is constant (in the absence of scattering interactions), ii) The size of the volume 
element, being associated with one of Poincare's integral invariants (under contact transformations), 
is also a constant of motion. Then, Liouville's theorem follows: Iff there are only conservative forces, 
then the phase space density is a constant of motion. 

Beginning with this theorem, we derive the evolution equation for / in terms of the relativistic 
proper time r as follows: 

= ^T/foP) = ^£S 4 Mt)-*]<5 4 [k(t)-p] , (2) 

i 

where we used the constancy of the density which is represented in terms of the particle (index i) 
trajectories {xi(r),pi(T)} in the first line, carried out the differentiation in the second, and employed 
the definition of the four- velocity and the equation of motion in the last, respectively; here ^(x) 
denotes the external or selfconsistent internal four-force(s). This is the relativistic Vlasov equation. 
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More generally, allowing for scattering of particles into and out of phase space volume elements, 
i.e. (d/dr)f(x,p) ^ 0, the Vlasov equation is replaced by the Boltzmann equation: 

-p-d x +T(x)-d p }f(x,p)=C[f}(x,p) , (4) 

where C[f] denotes the collision term. It can be derived from an analysis of the equations including 
the two-body densities and generally turns out to be a nonlinear functional of the one-body densities. 

A popular set of assumptions for this derivation is the following: i) Only binary or two-body 
collisions contribute (in a sufficiently dilute system); ii) Boltzmann's "Stosszahlansatz" according to 
which the number of collisions at x is proportional to f(x,p)f(x,p'); iii) the distribution function 
varies slowly on the scale of the mean free path, X m f p \V log f(x,p)\ <C 1. 

We remark that while the Vlasov equation is appropriate for systems with conservative forces 
only and, thus, describes nondissipative phenomena, the Boltzmann equation incorporates dissipative 
scattering processes, which lead to entropy production. 

In order to illustrate this, we introduce the simple relaxation time collision term , which could 
be derived rigorously as an approximation of a two-body scattering term but can also be seen as a 
phenomenological ansatz taking the dissipation into account: 

C\f\ = --(f-fo) , (5) 

where r r denotes the relaxation time parameter and /o the equilibrium one-body density, towards 
which the system will relax. 

A particular equilibrium solution of the Boltzmann equation with the above collision term can be 
obtained from the Jiittner distribution , 

fo(x,p) = eM-P(U-p + fi)} , (6) 

where the local parameters {(3 = 1/T, U^, denote the inverse temperature, flow four-velocity, 

and chemical potential, respectively. Instead of the exponential 'Boltzmann factor' one may also use 
a Fermi-Dirac or Bose-Einstein distribution, depending on the nature of the considered particles. 

In the example of a plasma of charged particles, including a static homogeneous neutralizing back- 
ground, the Lorentz force is J-^{x) = F^ v j v = (e/m)F fl ' / p u , in terms of the field strength tensor F^ v 
of external and/or selfconsistently generated internal electromagnetic fields and the particle electric 
current j^. Inserting this together with the Jiittner distribution, i.e. f = fo, into the Boltzmann 
Eq. (H) together with Eq. (||), we obtain the equation: 

l) [P{U-p + v)]+e(3F^U v = , (7) 



which constrains the parameters of Jiittner distribution. Thus, the simplest solution indeed is the 
global equilibrium distribution with a constant temperature, a global rest frame, and where the gra- 
dient of the chemical potential compensates the electric field F^, — c^/i = eF^ . 

2.2 Relation to relativistic hydrodynamics 

Turning to the observables to be described by the phase space density / introduced in Eq. (|l|), we 
define the particle (mass, charge) four-current, 



N»{x) = J d A pp»f(x,p) , 



(8) 
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and the energy-momentum tensor, 

T^(x)= I d*pp^f(x,p) , (9) 



as the first and second moment of the distribution function, respectively; it is important to recall here 
that we assume / to implicitly contain the usual on-shell constraint, i.e. the factors Q(p°)5(p 2 — to 2 ), 
as we discussed. 

These are the main quantities of interest in the hydrodynamic description of matter, where one 
integrates out the momentum space information contained in / fl6| 1. 

Indeed, it is straightforward to show that and T^ v obey the appropriate continuity equations 
related to mass (or charge) and four-momentum conservation. To begin with, using the Boltzmann 
Eq. (|j) together with Eq. (||), we obtain: 

fyJV" = Jd 4 pp- d x f = ~mj d 4 p (V • d p f + f —^j = -™ J d 4 p (/ - /„) = -dpSN" , (10) 

where we find a dissipative contribution 5N^ on the right-hand side, which vanishes only if the ordinary 
density equals the equilibrium density determined by /o- Similarly, one obtains: 

d^ = -m J d 4 pp v ^-d p f+ f -^j=m^{Jd 4 pf)- r ^{N^-Ni;) , (11) 

where the second term on the right-hand side is related to a dissipative contribution —d p 5T^ u to the 
energy-momentum tensor, while the first term presents the external or selfconsistent force density 
acting on the system; Nq is defined like N^, however, with / replaced by fg. 

More generally, based on the Chapman-Enskog method, standard forms of the dissipative terms can 
be constructed incorporating the transport coefficients of shear and bulk viscosity, heat conductivity, 
particle production/annihilation, and diffusion [16, fH|. 



Finally, we consider the entropy. It is an important quantity not only in equilibrium thermo- 
dynamics, but can be used, for example, to characterize the bulk properties of matter produced in 
high-energy collisions. In particular, it can be related more or less directly to the observed particle 
multiplicities [||, @, 16|- Here we define the entropy four-current: 



S^ X ) = -J d 4 pp»f(x,p)ln[f(x,p)/f (x,p)] . (12) 

Calculating as before, we obtain the entropy production formula: 

d x ■ S = - [ d 4 p (ln[/// ] + 1)(/ - /o) (13) 
r r J 

Since (lnx + l)(x — 1) > 0, for all x > 0, we recover Boltzmann's "H-theorem" : 

d x -S>0 , (14) 

expressing a positive entropy production which vanishes only in equilibrium, when f = f Q . 

This completes our overview of classical relativistic transport theory and some of its ramifications 
and we turn to quantum mechanics next. 

3 Quantum mechanical distribution functions 
- the density matrix and the Wigner function 

In order to motivate the necessity for a density matrix formulation in quantum mechanics, we recall 



Feynman's famous division [[ER 19], 



Universe = System + IZest , (15) 
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which says that doing physics and paying attention to the system requires precisely to separate off 
the rest of the Universe, i.e. the 'environment' of the system in modern terminology. 

Following the presentation in Ref. [18], we introduce generic coordinates x and y for the description 
of S and 1Z, respectively. Assuming the existence of corresponding complete sets of states or wave 
functions, {(f>i(x)} and {ipj(y)}, the most general state vector and normalized wave function of IA can 
be expanded, respectively, as: 

I*) =Z^jl^)hM . *(*»!/) = J2 c ij(y\ ( t>i}( x \i J 3} = ^2(h(v)<Ih(x) > ( 16 ) 

with complex expansion coefficients Cjj and functions Cj(y). 

Considering an operator A which acts on S only, it is given in terms of its matrix elements by: 

A = ^^|^)|^)(^l(0i'l , (17) 

ii'j 

i.e., it acts as a projector on 1Z. Following the rules, we obtain the expectation value of A in the state 
|*> of Eq.©: 

(A) = = ]T And rO/ = E A ViPw = TvAp = Tr pA . (18) 

i'ij' i'i 

It follows from its definition here that the density operator p, with the density matrix elements 
identified as pa' = {4>i\p\4>ii) , is a hermitean operator. 

Therefore, we may introduce a complete orthonormal basis {|i}} diagonalizing p: p = J2i w i\i) {i\ > 
with real eigenvalues Wi. Furthermore, choosing A = lg^ln, i.e. the identity operator, one finds with 
the help of Eq. (18) the sum rule: 1 = (A) = Tr p = J2i w i ■ Finally, choosing instead A = <8> 1-r., 

a similar calculation yields: Wi> = (A) = (^l^]*) = J2j li^jlii'ltyl 2 > . 

Abstracting from the present example, it is postulated that any quantum mechanical system is to 
be described by a hermitean density operator, 

p = ^2wi\i)(i\ , (19) 

i 

where forms a complete orthonormal set, the real expansion coefficients are non-negative, Wi > 0, 
and fullfill a normalization condition, J2i w i = 1- ~~ We remark that the density matrix was first 
introduced by von Neumann in 1932 in his by now famous book (2^[ . - Observables in particular and 
expectation values of operators in general are always to be calculated by: 

(i) = Tr pA = J2wi(i\A\i) . (20) 

i 

Consequently, the coefficients Wi are interpreted as describing the probability to find the system in 
the state 

As a matter of nomenclature, one distinguishes pure and mixed states of a system, the former 
being defined by a density operator which is a projector, p = p 2 <^=^> (wi* = 1, all other wi = 0), and 
the latter comprising all other cases. - We recall the case of a system in contact with a heat bath, 
which leads to the mixed state density operator of thermal equilibrium: 

p(P) = Z-\f3)J2e-P E "\E n ){E n \ , Z{[3) = Y. e ~ PEn > ( 21 ) 

n n 

where {|-E n )} denotes the complete set of normalized energy eigenstates of the system and Z is its 
partition function necessary to normalize the exponential 'Boltzmann factor'. All canonical thermo- 
dynamical relations can be easily derived from this density operator, see, for example, [18]. 
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Finally, we may address also in the present context the question how the system, i.e. its density 
operator, evolves in time. This is easily answered by expanding the eigenstates {\i)} of p in terms of 
the energy eigenstates, 

\i(0)) = \i) = J2\E n ){E n \i(0)) , (22) 

n 

which implies according to the Schrodinger equation: \i(t)) = exp(— iHt)\i(0)), given the Hamiltonian 
H. Then, the evolution of p follows: 

p{t) = e- lAt p(0)e +iAt . (23) 

Equivalently, we obtain: 

±m=wi&\ > (24) 

which differs from the usual operator evolution equation in quantum mechanics by an additional minus 
sign on the right-hand side. Both signs, of course, are consistent, since we have: 

(i) (t) = Tr p{t)A = Tr p^e+^ie"^ = (A(t)) , (25) 

using the cyclicity of the trace. 

The equations (^4|) and (^) present the problem of quantum transport theory in its most compact 
form. In particular, Sections 4-6 are devoted to explorations of various more detailed forms of these 
abstract results and their applications. 

In order to connect the density matrix formalism with the classical transport theory, we turn to 
the density matrix in coordinate and momentum representation. 

It helps to visualize the following by imagining a single particle to be described quantum mechan- 
ically. Then, in the coordinate representation, using corresponding single-particle wave functions, we 
have the density matrix elements: 

p(x',x) =^2wi{x'\i)(i\x) = ^2wi(pi(x')(p*{x) , (26) 



see Eq. (19). Then, we calculate immediately the probability to find the particle at x: 

P(x)=p(x,x) = J2w i \cl> i (x)\ 2 , (27) 

i 

in agreement with the standard rules of quantum mechanics. Similarly, in momentum representation, 

P(P,P) =X>i(p'N>(*b) = X>i&(pO#(p) , (28) 

i i 

which yields: 

Pip) = P(P,P) =X>i|&(p)| a , (29) 

i 

i.e. the probability to find the particle with momentum p. 

Now, considering a function O - representing some observable, for example - which is defined over 
phase space, we obtain its average value by calculating: 

= 1 dxdpO(x,p)f(x,p) , (30) 

involving the one-body probability density function /. The question then arises, whether there exists 
a corresponding quantum mechanical density function which yields the expectation value of operators 
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in the form of phase space integrals, generalizing Eq. (pQ). One possible answer is provided by the 



Wigner function pa, 19, 21 



W(x,p) ^Jdy e^' n p{x + y/2, x - y/2) , (31) 

i.e. a particular Fourier transform of the density matrix in coordinate space. We remark that the 
momentum appearing as an argument of W is the conjugate variable to the relative coordinate sep- 
arating the pair of wave functions which enter p, cf. Eq. (p6[); we also indicate here explicitly the 
^-dependence arising in the phase factor. 

Furthermore, even though the Wigner function presents a real distribution, due to W(x,p) = 
W*(x,p), it may oscillate and indeed does so in most cases. Nevertheless, we easily obtain the 
following results which match their classical counterparts: 

<j/ ' W(x,p) = p(x,x)=P{x) , (32) 



2ttTi 

dxW(x,p) = p(p,p) = P(p) • (33) 



This implies for the expectation values of functions of operators: 



(0(x)) = TvpO(x) = J^O(x)W(x,p) , (34) 

<0(p)> = TrpO(p) = J^O(p)W(x,p) , (35) 

which should be compared to Eq. (]30|). 

However, besides the fact that the Wigner function generally is not positive definite, the expected 
limitation of the classical/quantum correspondence of various (probability) density functions also 
shows up, when one considers operators of the form 0(x,p). In this case, operator ordering becomes 
an issue and the above formulae have to be generalized with care. On the other hand, as we shall see 
shortly in Section 4, the appropriate semiclassical expansion of the Wigner function transport equation 
does yield the transport equation for the classical probability density function, which we discussed in 
Section 2. 



4 Transport theory for quantum fields 

As an example of a field theory we choose the 0(4) linear a -model, which has a long history of 
applications in various phenomenological contexts. Recently it has been argued by Wilczek that it 
represents the QCD chiral order parameter for rif = 2 massless quark flavors Furthermore, 
it has been demonstrated by nonperturbative calculation that this model possesses a second order 
finite temperature phase transition between the spontaneously broken and symmetry restored phases 
1 23]. Most interestingly, the effective mass and coupling go to zero at the phase transition with well- 
determined critical exponents. This may influence the hydrodynamic behavior of matter described by 
this model in interesting ways, which we presently study. 

Here, our aim is to derive the quantum transport equations in the Hartree approximation. In 
particular, we will illustrate how an appropriate Wigner function can be introduced and how field 
theory aspects are related to our earlier considerations of classical transport theory. 

To begin with, the 0(4)-invariant c-model action is defined by: 

S$ - / d 4 * - - , (36) 
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with 4> — {<Pii ^3; <Pa) an d where the mass parameter is chosen with the 'wrong' sign, i.e. u 2 < 0. 
The corresponding potential, V{4>) = \l^4> 2 + A ( 2 ) 2 , is of the so-called 'Mexican hat' form which 
leads to spontaneous symmetry breaking at the tree level. 

We parametrize the four-component vector <j> in terms of three-component 'pion' and one-component 
'sigma' fields, 4> = & t each space-time point. Shifting the a'-field by the vacuum expectation 

value o"o of </>, which is determined by the minimum of the potential, 



dV(cf>) _ . 2 , 1 



(M 2 + ^A^=0 \a \ = J-6^/X, (37) 

we define the fluctuation field a = a' — oq. 

The effect of this parametrization combined with the shift of one component by the vacuum 
expectation value is easily seen in the corresponding Heisenberg operator equations of motion, 

<9 2 7? + m 2 n TT + ^vf 2 7f = , (38) 

3 2 a + m 2 a a + ^a a 2 + ^ycr 3 = , (39) 

which are obtained by varying the action S[<f)] with respect to </>, introducing parametrization and 
shift, and considering the fields as quantum field operators. Here we introduced the effective masses: 

ml = ^(a 2 + 2a a) , m 2 = ^(vf 2 + 2a 2 ) . (40) 

. We observe that for a — ► we have m 2 — > 0, while m 2 . —* Ao" 2 /3, yielding one massive 'radial' mode 
together with three massless 'Goldstone modes', in accordance with the Goldstone theorem. 

Taking the expectation value Eqs. (|38|) and (39), we observe that one-point functions, such as the 



mean field a{x) = (a(x)), generally are coupled to two-point functions, such as lim x ' ^ x {a{x)a{x')) , 
and higher n-point functions. The coincidence limit of these Wightman functions produces divergences 
necessitating a renormalization procedure, which we will discuss at the end of this section. 

In order to solve the resulting equations, one has to specify the density operator which enters the 
expectation values, e.g. (a) = Tr pa, omitting the operator signs used previously, cf. Section 3. This 
can be done, for example, on a fixed time-like hypersurface, which will be demonstrated for fermions 
in Section 6. Furthermore, we need to derive equations of motion for the higher n-point functions, 
since taking expectation values of nonlinear operator equations automatically generates an infinite 
hierarchy of equations, similarly as the Schwinger-Dyson equations for propagators or the BBGKY 
hierarchy in classical transport theory. 

The simplest nonperturbative truncation of the hierarchy of Wightman function equations is pro- 
duced by the Hartree approximation. It consists in factorizing the n-point functions into products of 
one- or two-point functions, properly taking into account all possible factorizations. For example, we 
obtain: 

(<7l<7 2 <7 3 ) = Ol^Ss) + ^2(^10:3) + ^3(0:10:2) + 010-203 , (41) 

where the subscripts refer to different space-time points. Note that the expectation value of any odd 
power of the proper quantum field vanishes, since we define a = a — a. - This approximation is known 
to be equivalent to summing all iterated bubbles ('superdaisies') in the Feynman diagram calculation 
of the vacuum effective action in a </> 4 -model [ 24 1 . 



Here we will furthermore assume that cross terms vanish, e.g. (vrc) = 0. It turns out to be 
consistent with this assumption to set (7?) = 0, since the classical pion field obtains no source term. 



In distinction, the Eq. (39) yields the Klein-Gordon type mean field equation: 



'd 2 + ml+ ^(V + a a + (a 2 )) a + ^a (a 2 ) = (d 2 + m 2 + 5m 2 ) a + J = , (42) 
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where the two-point function adds to the nonlinear classical force term shifting the effective cr-mass, 
which contains (vf 2 ), and to a source term, 5m 2 and J, respectively. For a homogeneous system, 
for 8m 2 ^> m 2 , and for sufficiently large (a 2 ), the mean field is (a) = —(Jo, indicating symmetry 
restoration with (a') = 0. 

Multiplying Eq. (^) by one more power of the field operator from the left, before applying the 
Hartree approximation as before, we obtain the two-point function equation: 



d 2 2 + ml 2 + Xa a 2 + -(a 2 2 + <a 2 2 }) ) (a lS2 ) = , (43) 



A 
2' 

where subscripts "1,2" refer to two different space-time points. Here we also used Eq. fl42| ) at point 
"2", multiplied by a\, in order to simplify Eq. (43) considerably. A similar equation follows from 
Eq. (^) for the 'pion' modes. 

The next step consists in defining a suitable Wigner operator, 

W ab (x,p) = J Aye-^y^ix + \y)Mx ~ \v) , (44) 

where $ = (7?, a) and $ = <i> — <I>. This should be compared to the quantum mechanical Wigner 
function introduced previously, Eq. (^|). We remark that it might be useful for some applications 
not involving vacuum properties to normal-order the field operators in the definition of W, see, for 
example, Section 6. 

Illustrating the usefulness of W, we write down the energy-momentum tensor for the 0(4)-model: 
T,„ = ((^)-M)-^[i(5$) 2 -^ 2 $ 2 -iA($ 2 ) 2 ]^ (45) 

+9»uJ d A p ^ 2 (W aa (x,p)) + ^\J dV (W aa (x,p)W bb (x,p'))^j , (46) 

to which must be added the purely classical terms plus mean field dependent interaction terms (oc A) 
involving W, which can all be further evaluated in Hartree approximation. 

Heading for the transport theory, we will express Eq. ( f43[) in terms of the Wigner operator. We 
introduce the abbreviation: 

M 2 (x) = ml(x) + \a a(x) + ^(a 2 (x) + (a 2 (x))) . (47) 

The two-point functions contained here, (ir 2 ) in m 2 and (cr 2 ), respectively, can be rewritten using the 
Wigner operator, e.g.: (c? 2 (x)} = / d 4 p(W aa (x,p)). Then we obtain instead of Eq. (|43|): 

1 i 
-d 2 -p 2 + ip-d x + exp(-- 



-d 2 -p 2 + ip-8 x + exp(--d x ■ d p )M 2 {x) ) {W aa (x,p}) = , (48) 



where the x-derivative in the exponential acts only on A4 2 . In the derivation of this result, and similarly 
for the 'pions', one frequently makes use of suitable partial integrations under the momentum integral 
from W, as well as expressing the shifted argument of A4 2 by a translation operator giving rise to the 
exponential 25]. 

The final step consists in adding to/subtracting from the complex Eq. (48) its adjoint. This yields 
the transport equation: 

p . dx _ I . d p )M 2 (x)) (W aa (x,p)) = , (49) 
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together with a generalized mass-shell constraint : 

4 ~ a 2 



—d 2 - cos(-d x ■ d p )M 2 {x) (W ao (x,p)) = , (50) 



where the appropriate powers of h are reinserted. We stress that the covariant transport equation 
alone does not suffice to determine the dynamics. 

A systematic semiclassical expansion of Eqs. (flS|)-(5C) in powers of h becomes feasible now. To 
leading order we obtain: 

(v-d x -^d x -d p M 2 {x)^{W (Ta {x,p)) = 0, (51) 

(p 2 -M 2 {x)) (W aa (x,p)) = , (52) 

which are indeed of the form of a classical Vlasov and mass-shell equation, respectively, discussed 
in Section 2.1. Note, however, that the effective mass or effective potential, appearing here contains 



contributions from the mean field, which have to be determined selfconsistently from Eq. ([42|), and 
'self energy' terms, which we will discuss shortly. 

We remark that effects of the higher-order ^-corrections to Eq. (|5l"l) and ([52] ) are largely unexplored, 
since numerically the corresponding higher derivatives lead to strong instabilities. This provides part 
of the motivation to develop an analytical approach solving the exact quantum transport equations, 
the first step of which is described in Section 6. 

Furthermore, we observe that the equations derived here do not yield any collision terms. This 
was to be expected, since the Hartree approximation, which is essentially equivalent to the Gaussian 
approximation in the Schrodinger functional approach, leads to an effective 'quantum' Hamiltonian 



evolution which is nondissipative 11]. Only an improved treatment of correlation terms, such as the 
last term in Eq. (^), will go beyond the collisionless Vlasov dynamics obtained here, see, for example, 
Refs.i|. 

Finally, we turn to the 'selfenergy' terms mentioned above and to the divergences caused by them, 
in particular. These terms involve integrals of the kind J d 4 p(W aa (x,p)). Let us consider the simplest 
case of a density operator projecting on the vacuum state, p vac = |0)(0|. Then, for a generic free scalar 
field, the vacuum Wigner function is: 

W vac {p) = Tr Pvac W(x,p) = (0\W(x,p)\o) = (2tt)- 3 <% 2 - m 2 ) , (53) 

where m denotes its physical mass, and the calculation proceeds, for example, by expanding the field 
operators appearing in the definition of the Wigner operator in terms of creation and annihilation 
operators. In this case, we obtain: 

d 4 p W vac (p) = J = I m , (54) 

with Up = {p 2 + m 2 ) 1 / 2 . This 'selfenergy' integral is quadratically divergent. 

Since divergences arise from field operators at coinciding points and, as a first approximation, 
independently of which point it is, we rewrite Eq. ([50]) accordingly, using Eqs. fllOD and (47): 



= (-p 2 + m 2 bare + -I m + -J d A k {W{x, k) - W vac (k)})(W(x,p) - W vac (p) + W vac (p)) (55) 
-p 2 + m 2 + ^J d A k W r (x, k)j (W r (x,p) + W vac (p)) , (56) 
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with W(x,p) = (W ua (x,p)) and the Lagrangian 'bare' mass rn 2 are = ^cr^, cf. Eq. (|40|). Here we left 
out all 'pion' and mean field contributions to keep things transparent. These can be added following 
the same strategy, while the gradient corrections would force us to keep W(x,p) unrenormalized at 
first sight. 

Note that we added and subtracted suitably the divergent term in Eq. 
Eq. (^) the renormalized mass and Wigner function, m r and W r , respectively. We complete the 
renormalization by identifying the physical mass: 

A f 

m 2 = m 2 + - d*k W r {x, k) . (57) 

With this, the Eq. ( |56| ) assumes a simple form, (p 2 — m 2 )W r {x,p) = 0, indicating that W r is (in our 
first approximation) an on-shell distribution: 

W r (x,p) = 5(p 2 - m 2 )f(x,p) , (58) 

such as a Bose-Einstein distribution with (weakly) re-dependent parameters, cf. Section 2.1. Neglecting 
this dependence altogether and assuming a simple thermal distribution turns Eq. ( |57| ) into a mass-gap 
equation with m 2 as input parameter and a resulting temperature dependent mass m 2 (T). 

A mass renormalization as performed here is sufficient to render the mean field Eq. ( |42"|) and 
the operator of two-point function Eq. (^) finite. The gradient corrections in the transport and 
generalized mass-shell constraint, however, indicate that the renormalized Wigner (or equivalently the 
two-point) function has to be be renormalized separately, subtracting as above, and thus closing the 
set of equations. Generally, the situation gets more complicated, when the x-dependence of W r has 
to be taken serious and, thus, the physical mass becomes space-time dependent, e.g. in the case of a 
system which evolves from strongly inhomogeneous initial conditions. We will not address these issues 
here, but turn to a related physical effect in the next section. 



55), thus introducing in 



5 Particle production by classical fields 

The color string, rope, or flux tube models are widely used in phenomenological descriptions of particle 
production, in high-energy nuclear collisions in particular The basis here is Schwinger's nonpertur- 
bative calculation of vacuum decay due to charged particle production in constant and homogeneous 



external electric fields [26, 27]. Considering the inhomogeneous and rapid evolution of the system 
during a heavy-ion collison, the basic assumptions of this picture are questionable. 

Since an appropriate generalization of Schwinger's result is still not available, one may as well 
consider the perturbative evaluation of particle production, however, for arbitrarily varying fields. 
This was originally discussed in Ref. pg] for QCD in abelian dominance approximation. There, also 
the ensuing modification of semiclassical transport equations, cf. Section 4, and the related vacuum 
polarization current were obtained. Here we briefly recall part of this calculation which is based on 
an 0(g 2 ) solution of the quark transport equation in external color fields. 

Simplifying the notation, we consider electrically charged Dirac fermions (mass m) in arbitrary 
electromagnetic fields. We define the Fourier transformed vacuum Wigner function: 

fvac(q,p) = |d 4 xd 4 ye^e-^(0|^(x + ^)V(x-^)|0) (59) 
= -(2ir) 5 5\q)5(p 2 - m 2 )e(-p°)(~f ■ p + m) , (60) 



which is a 4x4 spinor matrix. The final form here follows similarly as in the scalar case, Eq. (|53|), 
using a standard expansion of the field operators [|27[ . In distinction to the fermion Wigner function 
introduced in the following Section 6, Eq. (|69|), where the notation is explained in more detail, we 
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presently do not normal-order the field operators, since we are interested to study the response of the 
vacuum, i.e. the modification of f vac , due to external fields. 

Our task is to somehow solve the full quantum transport equation which determines the fermionic 
Wigner function. For the present study the transport equation based on the linear form of the Dirac 
equation is most useful || |25|]. - The quadratic form yields transport equations which are closer to 
the usual form, such as in the scalar field case studied in Section 4. - Presently, following Fourier 
transformation, we have: 

f i \ {' ' 1 

(l • {P + g?) - m J filiP) = 9 J ^-|4 7 • A(q')f(q -q',p- -q') , (61) 

where denotes the vector potential of the external field. This equation has to be solved together 
with the constraint : 

fHq,p) = i°f(-q,ph° , (62) 



which follows from the definition of the fermion Wigner function; equivalently, one may solve Eq. ( pT| ) 
simultaneously with its adjoint. 

The perturbative solution in powers of the coupling constant g of Eq. (|6l]) can be found interatively, 
starting with the zeroth order solution f vac as input on the right-hand side. 

Thus, at first order one finds the induced vacuum current . J^\(q) = J d 4 pTi'j fl f^(q,p), which 
can be written explicitly involving the first order (one-loop) QED vacuum polarization tensor and 
reproduces the known result [27, 28]. Furthermore, since there is a gauge invariant version of Eq. (|6l|) 



based on a modified definition of the Wigner function |g, |25|, 28 1, this type of calculation might allow 
to calculate the polarization tensor in a manifestly gauge invariant way. 

For our present purposes, the second order solution of the equations is most interesting. The main 



algebraic complication arises from Eq. (62), which has to be satisfied order by order. We do not give 



the lengthy expressions here, but turn to the calculation of the spectrum of the produced particles. We 
have to use reduction formulae, in order to relate the 'in' field operators defining the Wigner function 
to asymptotic on-shell 'out' particle states In this way, we obtain: 

2 ^ P ~¥p~ ^ ^(P) = 2 (2ir) 3 u r (p)(l-p- m)/ (2 )(0,p)(7 -p - m)u r (p) (63) 

= |^/d 4 g% -p>((g-p) 2 -m 2 ) 

• (%) • E(-q) - B(q) ■ B{-q) - \A{q) ■ (q - 2p)\ 2 ) , (64) 



where the electromagnetic fields and the vector potential enter; here uj p = ^p 1 + m 2 and u r (u r ) 
denotes the single-particle (adjoint) spinor wave function, respectively, in the normalization of Ref. |27f| . 
Also the last term of this result is gauge invariant, as can be demonstrated using the constraints of 
the integral and the fact that the particles are on-shell. Interestingly, with the external fields left 
completely general, still there is no spatial dependence of this spectrum. 

Finally, we calculate the vacuum decay rate 1Z per unit four-volume, using a more suitable inter- 
mediate form of the spectrum result, Eq. (164^ 

K = 2 /|^(2)W ( 65 ) 

= ^ l^^d\5\q-p-p')^q 2 g, u +p,pl+p^A^q)A^-q) (66) 

= ^/dV(g 2 -4m 2 ) (l-^H ( H-^-H (lEUDf-Wrtf) . (07) 
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which confirms the result obtained by other methods in Ref. [27|. Clearly, particle production is an 
electric field effect. 



For completeness, let us quote Schwinger's nonperturbative result [26|: 




^=^rE 4exp U^L , ( 68 ) 



which forms the basis to date of the phenomenological models mentioned in the beginning of this 
section. For reasons which we discussed, a future calculation bridging the gap between our perturbative 
calculation for arbitrary fields and the nonperturbative result for constant electric fields would be 
extremely useful. 



6 Fluid dynamics of relativistic quantum dust 

In this section, we will study the relation between relativistic hydrodynamics and the full quantum 
evolution of a free matter field p9j ]. In particular, we try to answer how a free fermion field and 
its energy-momentum tensor will evolve, given arbitrary initial conditions and especially those of the 
Landau and Bjorken models. 

In the absence of interactions, decoherence or thermalization may be present in the initial state, 
corresponding to an impure density matrix, but is followed by unitary evolution. We consider this as a 
"quantum dust" model of the expansion of matter originating from a high energy density preparation 
phase, which the Landau and Bjorken models describe classically ||] . 

Our approach is independent of the nature of the field, as long as it obeys a standard wave 
equation. To be definite, we choose to work with Dirac fermions and comment about neutrinos 
later. We introduce the spinor Wigner function, i.e., a (4x4)-matrix depending on space-time and 
four-momentum coordinates: 

W a p(x;p) = J T^e-^: foix + y/2)Mx ~ v/2) :> , (69) 

where the expectation value refers to the (mixed) state of the system; without interactions, the vacuum 
plays only a passive role and, therefore, is eliminated by normal-ordering the field operators. Note 
that the normalization of the Wigner function is a matter of convention, and we have chosen the most 
convenient one for us here. 

All observables can be expressed in terms of the Wigner function here. In particular, the (unsym- 
metrized) energy-momentum tensor: 



(: T^(x) :) = i(: ^(x)^ d„ i>(x) :) =tr 7 „ / d*pp„W(x;p) , (70) 
where d= \{d — d) and with a trace over spinor indices (conventions as in plj ) . Furthermore, the 



dynamics of W reduces to the usual phase space description, as in Section 2, in the classical limit [25]. 

The following study is based on the simple fact that the propagation of the free fields entering 
in Eq. (|69|) from one time-like hypersurface to another is described by the Schwinger function. It is 
the solution of the homogeneous Dirac equation, [27 • d x — m]S(x,x') = 0, for the initial condition 
S(x, x' , x° = x'°) = —ij°5 3 (x — x'). Thus, i/j(x) = i j d 3 x'S(x, x') r y°i/j(x'), and similarly for the adjoint. 
An explicit form is: 

iS(x, x') = iS(x -x' = A) = (17 • d A + m) f j- (e~ ik +- A - e~ ik - A ) , (71) 



(2vr)32^, 



where k± = (±u>k, k) and u>k = (k 2 + m 



2U/2 
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Making use of Eqs.© and (0), we relate the Wigner function at different times, t = x", x 



o ^/o. 

d 4 k 



W(x;p) = J ^ e - lk - x S ± (pi)5 ± (p^ J dVe w J dp' A( P +h°W(x';p'h A(p^ ) , (72) 

where pj^ =p ± |, ^(g) = ±5(q 2 - m 2 ) (for q°/\q°\ = ±1), A(g) = 7 • g + m, and p'^ = (p'°,p). 

The Eq. ( [72] ) implies that the Wigner function obeys a generalized mass-shell constraint and a 
proper free-streaming transport equation: 

h 2 

[p 2 -m 2 --d 2 x }W(x;p) = , (73) 

p-d x W(x;p) = , (74) 

separately for each matrix element. The reinserted fi indicates the important quantum term in the 
equations, which otherwise have the familiar classical appearance. 

Thus Eq. ( |7"2| ) presents an integral solution of the microscopic transport equations for a given initial 
Wigner function. Furthermore, a semiclassical approximation of the Schwinger function may be used 
to generate an integral solution of the corresponding classical transport problem. 

Next, we decompose the Wigner function with respect to the standard basis of the Clifford algebra, 

^(T^S^i/, i.e., in terms of scalar, pseudoscalar, vector, axial 
vector, and antisymmetric tensor components p7]]. The functions, J- = jtr W , V = — -A tr ^y 5 W , 
V M = jtr j^W , = |tr 757^ W , and S^ v = |tr cr^W , which represent physical current densities, 
are real, due to = ^Wj |^5| . They individually obey Eqs. (|7^ ) and (ff^). 

We assume V = = .Au, i.e., we consider a spin saturated system for simplicity. - From this 
point on, a corresponding study of (approximately) massless Standard Model vl^r neutrinos differs, 
and is simpler, since V M = A^ , while all other densities vanish identically; see, e.g., Ref. |30|1 . 
Then, using the 'transport equation' which follows directly from the Dirac equation applied to W, 
[7 ' (p + ^d x ) — m]W(x;p) = 0, and decomposing it accordingly, the following additional relations 
among the remaining densities are obtained: 

V^x;p) = ^Tix-p) , (75) 
p £ 

S^{x-p) = ^{j?d»-p»cQT{x-p) . (76) 

Note that S^" is intrinsically by one order in fi smaller than the other two densities. 

We conclude that presently the dynamics of the system is represented completely by the scalar 
phase space density which obeys the same transport equations as the full W itself. Using Eqs. ( |70|) 
and ([75|), we obtain in particular: 

p 

which is symmetric and conserved, d^T^ v {x) = 0, on account of Eq. (|74|), Furthermore, this implies 
the 'equation of state': 

(: T 00 (x) :) - ]T(: T u :) = 4m f d 4 p F(x;p) , (78) 
i=i J 



: T^(x) :} = 4m / d 4 p ^^T(x;p) , (77) 
J V 



which relates energy density and pressure(s). However, applying Eq. (73), we find that this relationship 



evolves in a wavelike manner, driven by off-shell contributions to the evolving T: 

dl{: TMx) :} = 16m / d 4 p {p 2 - m 2 )^(x;p) . (79) 



16 



H.-Th. Elze 



This differs from classical hydrodynamics with a fixed functional form of the equation of state. 
Eqs. (|77|)-(f79|) hold independently of the initial state, of course, if it evolves without further inter- 
action. 

Making use of Eq.tfrD in Eq. ([T?]), we now calculate the energy-momentum tensor at any time 
in terms of the initial scalar density. Employing the decomposition of the Wigner function and 
commutation and trace relations for the 7 matrices, as well as Eqs. (|73|)-(|76|), we obtain: 

(: T-»(x) :) = 8mJ...J P -^(p 2 + ^ (pV + jf 2 ) - ^ ((AT - P ^ 



+P%' k 2 + ^(k )V)jHx'-,p, P ' ) , (80) 



where J ■ ■ ■ J = (2vr)~ 3 J d 4 p J d A k e~' lk ' x 5 ± (p'^)5 ± (p^) J d 3 x'e lkx ' J dp' ; we also made use of partial 
integrations and the ^-function constraints. The three terms on the right-hand side stem from the 
scalar, vector, and antisymmetric tensor components of the initial Wigner function, respectively. 

If the initial distribution is an isotropic function of the three-momentum, then T^ v is diagonal at 
all times, implying that the absence of flow in the initial state will be preserved. 

Indeed, we expect the (non-)flow features of the initial distribution to be preserved during the 
evolution, due to the absence of interactions. Kinetic energy from microscopic particle degrees of 
freedom will not be converted into collective motion. An interesting question is, how the classical hy- 
drodynamic acceleration of fluid cells due to pressure gradients arises in our present model after coarse 
graining pCj , O , [i"2| , [l3|| . This has not been studied yet. Recalling earlier work on the hydrodynamic 
representation of quantum mechanics, e.g. Refs. [Si]], and recently deduced classical fluid behavior of 
quantum fields in WKB approximation ||, however, we study the full quantum evolution here. 

We are particularly interested in the exact evolution of T^ u , assuming a particle-antiparticle sym- 
metric initial state. This is believed to hold, for example, close to midrapidity in the center-of-mass 
frame of central high-energy collisions ||, ||] . It implies that the initial T is an even function of the en- 
ergy variable, T{x ';p,p' ) = T{x';p, —p'°). While Eq. (|80|) allows general initial conditions, we follow 
the implicit on-shell assumption in classical hydrodynamic models: 

F(x';p,p'°) = (2vr)- 3 m 5{p' 2 - m 2 ) (e(p' )F{x';p,p'°) + Q(-p'°)F(x';p, -p' )) . (81) 



Fermion blackbody radiation is described by F(x';p,p'°) = f(p'°/T(x')), where T denotes the local 
temperature, and with f(s) = (e s + l) -1 ; this is easily illustrated with the help of Eqs. (|77[) and (|8T|). 
Implementing Eq. ([Si]), we obtain the simpler result: 

(: *~(.) =) = / ^TSni ^os[Mf-f<)] , 

+ uJ-) 2 — k ) cos[(w + — uj-)t] — — uJ-) 2 — k 2 ) cos[(w + + oj_)t]} , 

where t = x° — x'°, uj p = (p 2 + m 2 ) 1 / 2 , oj± = ((p ± k/2) 2 + m 2 ) 1 / 2 ; furthermore, p° = \\oj+ =b 
with "+" when multiplying the first and "— " when multiplying the second term of the difference, 
respectively. Depending on geometry and initial state, further integrations can be done analytically. 

Consider a (1 + l)-dimensional system for illustration, assuming that the particles are approx- 
imately massless, i.e. uj p ~ \p\, and that F is even in p (no flow). Specializing to a kind of 
Landau initial condition , the distribution is prepared on a fixed timelike hypersurface at t = 
Q. We find the ultrarelativistic equation of state for the only nonvanishing components of T^, 
e = T 00 = T 11 = P {d = 1 + 1), which are calculated as a momentum integral following Eq. (|32^): 

T 00 (x,t) = 2j^\p\(F(x-t;\p\)+F(x + t;\p\j) (83) 
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= I(t°V-M = 0)+T 00 (x + M = 0)) , (84) 

i.e., a superposition of wavelike propagating momentum contributions in accordance with Eq. (|79|). 

Similarly, a kind of Bjorken initial condition can be specified on a surface of constant proper time 
Q. A transformation of Eq. (|33|) to space-time rapidity and proper time coordinates yields: 

T 00 (y,r) = 2j ^[^(-roe"^^^; \p\) + F(r e^ 2+ln ^; \p\j) , (85) 

since x = rsinhy/2 and t = rcoshy/2 (r > To > 0). 

Our results for free fermions show the free-streaming behavior of classical dust , associated with 
the independent propagation and linear superposition of the momentum contributions to the scalar 
component F of the Wigner function here. In particular, the shape function of each mode is preserved 
and translated lightlike (with dispersion for massive particles). Due to the assumed momentum sym- 
metry, the initial distribution will separate into two components after a finite time, travelling into 
the forward and backward direction, respectively, with a corresponding dilution at the center. For 
example, an initial distribution of Gaussian shape will separate into two corresponding humps. 

We recall that T^ v being diagonal implies the absence of ideal hydrodynamic flow, given e = 
id — 1)P- This does not depend on whether the initial state is on- or off-shell, see Eq.(|80|). Therefore, 
any hydrodynamic behavior must be the effect of a pecularity of the semiclassical limit P], of coarse 



graining |Ll| , |l2j , 13], or of interactions [p^| , or a combination of these. 



Despite the apparently classical evolution, however, all initial state quantum effects are incorpo- 
rated and preserved. If the initial dimensionless distribution F has a dependence on products of 
momentum and space-time variables, which is characteristic for matter waves, such terms invoke a 
factor 1/h. Similarly, if it is thermal (T) but includes the finite size (L) shell effects or global con- 
straints, then there are quantum corrections involving LT/h (&b = c = 1) j33[|. For typical values 
of LT/h ~ 1 these latter corrections are known to lead to corrections on the order of 30% in the 
thermodynamical quantities. They have not been included in semiclassical transport or classical hy- 
drodynamic models of high-energy (nuclear) collisions, but may be large. Here the quantum dust 
model provides a valuable testing ground to assess the importance of these quantum effects. 

Let us summarize briefly the results and perspective of this section: 

• Based on the Schwinger function, we obtain the solution of the free quantum transport problem 
to quadratures for arbitrary on- or off-shell initial conditions. 

• In 1+1 dimensions free fermions, i.e. their observables embodied in the energy-momentum 
tensor, collectively behave like classical dust showing no flow effects, if there is no flow present 
in the initial condition. Corresponding analytical results in three dimensions can be obtained 
and will be discussed elsewhere. 

• The method presented here may lead to an efficient way of treating interacting particles. Espe- 
cially when a low-order perturbative expansion is meaningful, interactions could be incorporated 
in a multiple scattering expansion. With free propagation in between scattering events, a con- 
sistent and conceptually simple space-time description of transport phenomena seems feasible. 

Thus we conclude our introductory lectures on relativistic quantum transport theory, which should 
convey some of its basic concepts and hopefully will lead to some of the interesting topics for further 
study. 
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